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A proof that continuous time quantum walks are universal for quantu m computation, using un- 
weighted graphs of low degree, has recently been presented by IChilda [PRL 102 180501 (2009)]. 
We present a version based instead on the discrete time quantum walk. We show the discrete time 
quantum walk is able to implement the same universal gate set and thus both discrete and contin- 
uous time quantum walks are computational primitives. Additionally we give a set of components 
on which the discrete time quantum walk provides perfect state transfer. 
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t_i ' I. INTRODUCTION 



Quantum computers offer the promise of fundamentally faster processing based upon quantum mechanical prop- 
erties. Although a physical device of a useful size is still to be built, m any qu ant um algo rithms have already been 
discovered. The most important of these are the algorithms introduced bylShgJand[G roverl which can factor integers 
and search an unsorted database respectively, significantly faster than the best known classical algorithms [l|, H| . 
' Quantum walks were initially introduced in both continuous H and discrete time, in direct analogy with their 
classical counterparts, and have since been studied extensively Q. In the same way that classical random walks are 
used in computer science for algorithm design, many quantum algorithms have been developed based upon quantum 
walks, with varying speed ups over the best known classical algorithms for the same problem, Q. These solve the 
problems using two different approaches: hitting times and searching. In hitting time problems we start from a 
Q" 1 specific vertex and want to get to another as quickly as possible. These problems have yielded the largest speed up, 
including exponential speed ups over the classical case, 043 ■ Searching for an entry in an unsorted database is a 
| classically time-consuming problem taking on average a time of O(N) to search a set of N entries. Grover's algorithm, 
■ 0, improves on this to 0(\/~N) by using a technique known as amplitude amplification. The same speed up can be 
■^j" ' obtained using a q uantum walk method on various structures |10l 4l2j. 

In [H|, IChildsl extends the original results of iFevnmanl jl4| to show a co ntinuou s time quantum walk, on an 
unweighted graph of bounded degree, is universal for quantum computation. IChildsl gives an explicit construction 
that converts a standard gate model computation into a graph, on which a continuous time quantum walk executes 
an algorithm by traversing the graph. In this paper, we show the equivalent construction of a universal gate set using 
the discrete time quantum walk in place of the continuous. This confirms that both the continuous and the discrete 
time quantum walks can be regarded as computational primitives. The construction requires an exponentially large 
graph for the size of the input as we require 2™ wires for an n qubit input. The quantum walk takes place on this 
N- vertex graph just as the continuous time walk does in the construction by Childs (l3| . It is already known that a 
quantum walk on an N- vertex graph can be simulated efficiently by a universal quantum computer using poly (log N) 
gates, provided there is a simple rule to compute the neighbours of any vertex [7[. Thus, by performing the quantum 
walk on a quantum computer, the binary encoding brings the resources required back to the expected level. 

Our construction for the universal gate set in discrete time is similar to [H[ but has maximum degree, d, of eight 
at any vertex as opposed to three in the continuous case. The continuous time walk can easily be propagated in one 
direction with no reflection at the vertices. The discrete time walk is not so straightforward, it can only be propagated 
in one direction by using a specific coin corresponding to the a x operation. Using this coin restricts the graph to 
vertices of degree two, providing no way to construct higher degree structures. Thus we must use a double-edged 
wire to accomplish direc tional pr opagation. This solution has its roots in the connection between the continuous 
and discrete time walks. IStrauchl [l5| has shown that, as we take the continuous limit of the discr ete time walk on 
the line, we actually get two copies of the continuous time walk propagating in opposite directions. IChildsl [l6| later 
showed a di rect cor respondence between the discrete and continuous time quantum walks on arbitrary graphs. In the 
same work, IChildsl shows how a discrete time walk can be used, at its limit of small eigenvalues, to approximate the 
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continuous time walk. He uses this 'lazy' quantum walk approach to allow the discrete time walk to propagate in the 
same way as the continuous. This same approach could be used in this case to allow the computation to be performed 
on the same structures defined in [l3j]. However, this would require the discrete time walk to approach the limit at 
which it is doing very little at each timcstep. This would then increase the overhead required to allow completely 
deterministic computation. 

We begin by describing the discrete time quantum walk briefly in section HH then move on to show structures on 
which the discrete time quantum walk will allow perfect state transfer in section IIII1 These structures allow us to 
construct the elements we need to perform computation. Section [IVI shows the universal gate set we choose and how 
we implement these using the discrete time walk. Section [V] describes how we can link these gates and structures 
together to form any quantum circuit, and elaborates on how this is efficient, despite the size of the graph being 
exponential in the number of gates required. Finally, in section IVII we discuss our findings and the differences with 
the continuous time construction of [13j ■ 

II. DISCRETE TIME QUANTUM WALK 

Consider a classical random walk on a line in which a walker starts at a specific position and, depending on the 
outcome of a coin toss, moves either left or right. The outcome after many runs is a binomial distribution about the 
starting position with a spread (quantified by the standard deviation) of \fi where t is the number of timesteps. A 
discrete time quantum walk is the direct analog of the classical walk with the walker replaced by a quantum particle 
carrying a two state quantum system for the coin. The coin toss is effected by a unitary operator. Although this is 
now deterministic, if we were to measure the coin we would get a random output as in the classical case. We start 
the quantum walker at the origin and act upon it with a unitary operator for the coin toss, followed by a conditional 
shift operation (to obtain the movement of the walker) at each timestep. We write the basis states of the walker as 
an ordered pair, | x,c), denoting the position of the walker, x, and the state of the coin, i.e., heads (c = 1) or tails 
(c = 0). The simplest unitary operator is the Hadamard operator, H, which acts on the state as 

H \x,0) = -L(\x,0)+\x,l)) 

H\x,l) = -j=(\x,0)-\x,l)). (1) 



The shift operation, S, acts thus 



S\x,0) =| 1,0) 

S\x,l)=\x + l,l}- (2) 
The first three timesteps starting at the origin are as follows 

(SHf | 0,0) = (SH) 2 S^=(\ 0,0) + | 0,1)) 
= (Sff) 2 -^(|-l,0)+|l,l)) 
= (SH)S±(\ -1,0) +|-1,1) + | 1,0) -| 1,1)) 
= SH^(\ -2,0) + | 0,1) + | 0,0) - | 2,1)) 

= S^=(\ -2,0) + | -2,1) + | 0,0) - | 0,1) + | 0,0) + | 0,1) - | 2,0) + | 2,1)) 

= -L(|-3,0) + |-1,1) +2 | -1,0) -| 1,0) + | 3,1)). (3) 

As the walk progresses, quantum interference occurs whenever there is more than one possible path of t steps to the 
position. This can be both constructive and destructive, as shown in eq. which causes some probabilities to be 
amplified or decreased at each timestep. The walk on the line has been solved analytically [l7l . [l8j where it was first 
remarked that the quantum walk spreads quadratically faster than the classical one. 

The choice of operator at each vertex can greatly affect the dynamics of the walk and its propagation across the 
structure. A bias can be introduced 19], for d = 2, this is done using a generalisation of the Hadamard operator, 



FIG. 1: A general graph, G, with 5 vertices and 6 edges. It is undirected and has vertices of varying degree. 
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FIG. 2: An example of a portion of the walker reflecting back upon itself. This is a single degree 4 vertex with the Grover coin, 
eq. ((5j), operating on the incoming amplitude a. 
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where 8 is the bias in the coin. Setting this to 6 = h returns the standard Hadamard operator, eq. ([TJ. Similarly, 
the choice of the walker's initial state is also important, unlike in the classical random walk. A good review of these 
effects can be found in [20]. 

For universal computation we need a quantum walk on a more complex graph. In graph theory, a general graph, 
G, is an ordered pair consisting of a set of vertices, V, and a set of edges, E, which link the vertices. The number of 
edges incident on a vertex is the degree of the vertex, fig. [1] shows a small general graph which has vertices of varying 
degree. It is also undirected, meaning that the edges allow movement in both directions. A different operator (coin) 
is needed at vertices with d > 2 in order to act on the entire state space [2l|, [22| • 

The Grover coin can be extended to any degree at a vertex, 
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where d is the degree of the vertex and Id is the identity matrix of the same dimension. 

Coins of degree four will be needed at most of the vertices in our computational graphs. The Grover coin in four 
dimensions, eq. ([5]), reduces to 
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Using these higher dimensional coins can cause the walker to be reflected back upon itself with some probability. Figure 
[2] shows this reflection for a vertex of degree d ~ 4. In quantum walk search and other quantum walk algorithms 
this can be useful to provide interference. However, here we need to ensure the walker moves in one direction only, 
from left to right, so it most resembles the circuit model of quantum computation. We show how we accomplish this 
forward only propagation in section IIVI 



III. PERFECT STATE TRANSFER 



As a preliminary to our quantum computation scheme, we discuss structures on which perfect state transfer can 
be achieved using the discrete time quantum walk. Perfect state transfer has been investigated in the context of spin 




Perfect state transfer 
after 12 timesteps 



FIG. 3: Cycle of 8 vertices which gives perfect state transfer from the initial vertex to the opposite vertex after half of the 
period, 12 timesteps. The entire state returns to the initial vertex in a full period, 24 timesteps. 



chains by lLandahl et al.l 28L 29] . The propagation of the state through spin systems follows the same dynamics as a 
continuous time quantum walk. Perfect state transfer can occur on chains of length 2 or 3, hypercubes of any size, 
and chains with different coupling strengths engineered to optimise state transfer. 

The closely related properties of instantaneous mixing and periodic cycles have been studied in detail for q uantum 
walks. For the continuous time quantum walk, instantaneous mixing has been investigated by iTamon et al.l [24T - [27| . 
They showed [24| this is achieved on cycles of 2, 3 and 4 vert ices only. For the discrete time walk, slightly larger 
cycles show exact periodic behaviour. iTravaglione and Milburnl [23j j showe d that a cycle of 4 vertices has a periodicity 
of 8 timesteps after which the entire state returns to the starting position. iTregenna et al.l [2p| showed more periodic 
cycles exist, cycles of 2, 3, 4, 5, 6, 8 and 10 were shown numerically to be periodic by varying both the bias and phase 
in the coin. Perfect state transfer occurs at half the periodic cycle for even cycles, where we obtain the entire state 
at the opposite point of the cycle as shown in fig. [3] 

For our case of using the walk for computation, wc require the walk to travel perfectly in a single direction. On the 
structures mentioned, the quantum walk travels around the cycle in both directions and interferes to produce perfect 
state transfer. Using a completely biased coin, 
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we can make the state transfer perfectly around the cycle in a single direction. However, if we then try to attach 
another structure to the cycle, this periodicity is broken in both cases. The Grover coin, eq. ([5]), can be used to 
overcome part of this problem at vertices with an equal number of input edges as output edges. For any vertex of 
even degree it will transfer the entire state from all the input edges to all the output edges provided the inputs are 
all equal in both amplitude and phase. These results led us to the designs that work for universal computation. 

IV. UNIVERSAL GATE SET 

We now show how we construct a universal gate set with the discrete time quantum walk. Although the gate set 
we implement is the same as in [13j, the structures used to propagate the discrete time walk are different. The gate 
set used is the standard universal set comprising the controlled-not (C-NOT) gate, 



C 



NOT 



"1 








0" 





1 

















1 


.0 





1 






(8) 



the single qubit Hadamard, 
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and phase shift gates (we implement the specific phase shift known as the | gate), 



1 




(10) 



5 



rv n n n 

L-t u u u 

XIX 

a o oo 


Graver 


U L-fc u u 

XIX 

o a o o 


Shift 


n n n 

U U ' X u 

XIX 

o o a o 


A 




B 




c 



FIG. 4: Grover and shift operation acting on a vertex of degree d — 4. Section A shows the initial state, B shows the state 
after the Grover coin is applied and finally C is after the shift operation. 

FIG. 5: Basic wire used to propagate the quantum walk from left to right only. At a vertex of degree d — 4 the Grover diffusion 
coin is used. The initial state is split across the pair of edges in the wire, eq. 



These gates create a universal set that c an impl ement any quantum computation [301 ] . 

In order to represent quantum states, IChildsl defines his computational basis states as quantum wires. The other 
gates required for universality are then attached to wires and used to connect them together. The computation is 
represented as a quantum walk on these wires and structures, where the computation flows from input to output (left 
to right in our diagrams). Note that this encoding is not meant to be implemented directly. The wires represent 
computational basis states rather than qubits, thus the model does not represent a physical architecture. Instead, 
the underlying graph structure created would be used to help 'program' a quantum computer. We first show how to 
construct a simple wire along which the quantum walk will propagate naturally in one direction. We use two edges per 
wire to ensure no reflection occurs at a vertex. We distribute the walker across the two edges which then recombine 
at the next vertex. As the split is equal, the Grover coin in effect moves both halves to the output edges of the vertex. 
Figure 2] shows this operation and the 'shift' to the next vertex in explicit steps. The Grover diffusion coin, eq. ([5]), 
is used at each vertex of degree d = 4. The initial and final vertices are in effect degree four if we include other edges 
attached to either end. Figure [5] shows the basic wire we use. The computation would start with the amplitude at 
the initial vertex spread equally across the pair of edges in a wire. For example, the state a | 0) + f3 | 1), where 
| a | 2 + | ft | 2 = 1, would be split thus, 



I ^initial) = -j= [a | 0} a + a | 0> 6 + P | 1)„ + P | l) b ] , (11) 

where the subscript a refers to the top line of the wire and subscript b is the bottom line. The walk propagates left 
to right on the wire deterministically, in this case reaching the incoming edges of the final vertex in four timesteps. 
These wires form the basic connections in the computation. 

The simplest gate to construct is the C-NOT. It is trivial to implement by just exchanging the wires of the second 
qubit. The C-NOT gate is shown in fig. [5] and shows how the second qubit is flipped but the first qubit is untouched. 

The phase gate, eq. (|10p . requires the addition of a relative phase to one wire or computational basis state in relation 
to the other. To accomplish this, but still have only one coin operator for each vertex of the same degree, we modify 
the basic wire and add a phase factor, eX to it, 



Thus, as the walk propagates along a basic wire, it now picks up a phase of e 1 ^ each time it passes through a vertex 
of degree d = 4. For the wires shown in fig. [5] the walker would pick up a phase of e 5l< ^ as it reaches the final vertex. 
The phase added here is arbitrary and can be set to any value so long as it is set to the same value for all vertices of 
degree d = 4. As we are looking to implement a tt/8 gate we set it as follows: (f> = — 7r/4. In order to add a relative 
phase of 7r/4 between the | 0) and | 1) wires we insert the structure in fig. [7] into the graph at the required point. In 
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FIG. 6: Structure used to implement a C-NOT gate. In this case the first qubit is the control and the second is the target. 
The target qubit 's wires are interchanged and the control qubit is left unaltered. The dotted lines represent wires passing 
underneath the solid lines - there is no interaction between these wires. 
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FIG. 7: Phase gate structure. The d = 2 Grover coin, eq. (|13|l . is used at the vertices of degree d = 2. The | 1) wire will pick 
up a phase of e"* relative to the | 0) wire. In our construction we actually obtain the operation corresponding to a phase of 
e'~ as we set <f> = — 7r/4. 




FIG. 8: Hadamard gate structure. Sections A and C add a relative phase of i to the | 1) wire. The structure adds a global 
phase of 37r/4 to the wires. 



this structure there are also vertices of degree d = 2, at which we use the Grover coin at its limit of degree d = 2, 
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We add no phase to eq. (fl~3|) so as the walker passes through these vertices no phase is picked up. In fig. [Jj the walker 
propagates along the | 1) wire and picks up a phase of e~ i7r in four timesteps as it only passes through four vertices 
of degree four. However, the | 0) wire picks up a phase of e _5l T in the same number of timesteps as all its vertices 
are of degree four. Relative to the | 0) wire, the | 1) wire will pick up a phase of e 1 ^. Therefore, using the structure 
described here we obtain the operation in eq. (flU)) . 

The last gate in the universal set is the Hadamard gate. This requires an interaction between the two computational 
basis states. The structure we use to perform this operation is shown in fig. [8] This looks complex in relation to 
the other gates we have shown so we break it up to explain it more clearly. Sections A and C of the structure are 
each two phase gates giving a relative phase of i to the | 1) wire before and after the main section of the gate (B). 
Section B of the structure combines the two inputs from the | 0) and 11) wires and then splits this across the outputs 
equally. The structure here is similar to the basis changing gate in [13j. In order to obtain the desired operation on 
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this structure we have designed a coin for vertices of degree d = 8: 
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This operator combines the complex Hadamard operator, 
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in a tensor product form, 
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with the top two and bottom two rows of the Hi £g) Hi matrix rearranged. This rearrangement ensures the outputs 
come out in the same order as the input states. This gate adds a global phase of 37r/4. The phase gates at the start 
and end of the central section give us the relative phase of —1 required for the Hadamard operation. We note here 
that the choice of where to place the phases in this construction is arbitrary. The same result can be achieved by 
using the degree four Grover coin, eq. with no phase at vertices of degree four and the degree two Grover coin, 
eq. (fl"3|) . with a phase of ir/4 at vertices of degree two. However, by placing the phases on the Grover coin in the 
fashion we described previously, eqs. (|T2l [T3")) . means that the global phase added by the Hadamard gate, eq. (fl~4"|) . 
corresponds to the phase added by a wire of the same length. 



V. CONSTRUCTING QUANTUM CIRCUITS 



Thus far each gate we have described only acts on one or two qubits. However, non-trivial quantum computers 
involve many qubits. We now describe how to link these wires and structures together to form larger circuits. Figure [TU1 
shows the underlying graph structure of the circuit in fig. [UJ The graph structure is obtained by connecting together 
wires and structures so that the walk flows from left to right. For this reason, we designed our wire and structures 
with both input and output vertices of degree four, thus making it simple to link them together. The initial state of 
the computation is set on all or a subset of the vertices on the left hand side of the graph, with the amplitude at each 
vertex split across the incoming edges. This initial state can be thought of as the first column of vertices in the graph 
structure in superposition, with each subsequent column of vertices representing a further timestep. For example, in 
fig. [TU] this column of vertices is the set prior to the Hadamard structures. The walker is propagated across the graph 
structure, from left to right deterministically, for the required number of timesteps. We therefore do not require the 
addition of momentum filters or separators as in the continuous time case. Our structures all propagate the walker at 
the same speed, meaning output from the wires will be synchronised throughout the computation. Finally, the walker 
picks up a global phase of — tt/4 per vertex that is not part of a gate that changes the phase, so all the wires also 
stay synchronized in phase. Thus, we know with certainty that, after the required number of timesteps, the walker 
will have a distribution over just the output vertices on the right hand side of the graph. Once the computation has 
been completed, we measure the output vertices. We will find the walker at just one of these vertices, representing 
the output of the computation. 

The graph structure in fig. [TUJ is clearly larger in size than its equivalent representation in the circuit model, 
fig. [UJ In fact, for a general n-qubit computation the equivalent graph will have 2™ wires, one for each combination 
of computational basis states. Similarly, we require more gate structures than in the circuit model. Single qubit 
structures are repeated 2" _1 times and for the C-NOT gate we need 2 n ~ 2 structures. As an example, we can see the 
phase gate acting on qubit 3 in fig. [UJ is repeated 4 times in the underlying graph structure of fig. I1Q[ one for each 
combination of wires involving qubit 3. Although this seems as though we would lose any form of quantum speed up 
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FIG. 9: Quantum circuit on three qubits. A Hadamard operation is performed on qubit three followed by two C-NOT gates. 
Finally a phase gate is applied on qubit three. The underlying graph of this structure is shown in fig. [10] 
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FIG. 10: Graph to represent the quantum circuit in fig. [9] The Hadamard structure, H, is the same as in fig. [8] The dotted 
lines represent wires passing underneath the solid lines - there is no interaction between these wires. 



due to the exponential number of gates required in the underlying graph, this is not the case. Consider simulating a 
classical random walk on a classical computer, the iV-vertex graph is represented in log 2 N bits of memory with each 
vertex having a unique binary number as a label. In a similar fashion, if we simulate a quantum walk on a quantum 
computer, the N- vertex graph can be represented by log 2 N qubits. Therefore if we encode our graph using qubits, 
we can describe the 2 n wires in just n qubits. By manipulation of a single qubit we can affect all combinations of 
wires associated with that qubit. As the state moves across the graph, the adjacent vertices must be established. 
In complex graphs the description of the graph and its connections is often exponential in size and an oracle must 
be used to store it Q . The graphs produced here are of bounded degree and have a regularity stemming from the 
repetition of gate structures on combinations of wires involving a specific qubit. Due to the labelling of the wires, 
we know where to place each structure based on one bit in the label, thus we can efficiently describe the graph. For 
example, consider the second C-NOT gate in fig. [9j which operates on qubit 3 with qubit 2 as control. We can see 
from fig. [TU] that it is easy to identify where the C-NOT structures should be placed. The labelling of the wires shows 
that the middle bit determines which combinations of wires relate to the 2nd qubit having a value of 1. Similarly, we 
can also identify which wire it should link to by the last bit in the labelling scheme, it needs to be flipped relative to 
the original wire, i.e. | 010) links to | 011). 
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VI. DISCUSSION 

In this paper we have described an alternative to the construction in [l3[ using the discrete time quantum walk. This 
shows the discrete time quantum walk is universal, therefore any quantum algorithm can be reformulated as a discrete 
time quantum walk algorithm. It also confirms that the discrete and continuous time walks are both computational 
primitives and thus computationally equivalent. This equivalence is dependent on the number of steps in both cases 
to be of the same order. Our gate constructs require twice the number of edges compared to the continuous time 
case but the same number of wires. Our phase gate requires an additional timestep in relation to the continuous time 
phase gate construct. The number of timesteps required for a computation is also the same as the continuous time 
case but with a small overhead depending on the number of phase gates required. 

Another difference in the two constructions is the degree of the graphs produced. In the continuous time case the 
maximum degree of any vertex in the graph is three. In the discrete time case we use vertices of higher degree to 
ensure directional propagation. In most of the structures this is a doubling of the degree at a vertex, as shown in 
the case of the basic wire and the phase gate structure. The Hadamard structure we propose here howe ver, doe s not 
follow this doubling. It would seem reasonable, from the equivalent degree three structure in the work bv lChildsl . that 
it may be possible to decompose our Hadamard structure into one with degree six vertices. The doubling of degree 
at vertices would then correspond directly to the continuous case. 
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